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Abstract 

^ i We review different computation methods for the renormahsed energy momentum tensor of 

(N ■ a quantised scalar field in an Einstein Static Universe. For the extensively studied conformally 

coupled case we check their equivalence; for different couplings we discuss violation of different 
. energy conditions. In particular, there is a family of masses and couplings which violate the weak 

and strong energy conditions but do not lead to spacelike propagation. Amongst these cases is 
that of a minimally coupled massless scalar field with no potential. We also point out a particular 
coupling for which a massless scalar field has vanishing renormalised energy momentum tensor. 
We discuss the backreaction problem and in particular the possibility that this Casimir energy 
^ • could both source a short inflationary epoch and avoid the big bang singularity through a 

D , bounce. 

> '■ 

■ 1 Introduction 

H m . . 

I The Casimir effect [T] is a manifestation of the vacuum fluctuations of a quantum field. It was first 

considered in systems with boundaries, and it is known that the Casimir force is highly sensitive 
to the size, geometry and topology of such boundaries. In particular it may change from attractive 
to repulsive depending on such shape |2] ■ But the Casimir force is also present in systems with no 
boundaries and a compact topology, since the latter imposes periodicity conditions which resemble 
boundary conditions. 

If our universe is either open or flat, with non-trivial topology, or closed, every quantum field 
living on it should generate a Casimir type force, which led many authors to study the Casimir effect 
in FRW models (see |2] and references therein for a review). In the case of a spherical universe, most 
computations of the Casimir energy, or more generically, of the renormalised stress energy tensor, 
have focused on conformally coupled scalar fields. For instance, a massless conformally coupled 
scalar field, the electromagnetic field and the massless Dirac field on an Einstein static universe, 
have been considered in |B1 [7j; their Casimir energies have been shown to be of the form a/R^, 
with a, respectively, l/(4807r2), ll/(2407r2) and 17/(192077^). Note that all of these are positive. 
Since these are all conformally coupled fields, they have equation of state p = p/3, and obey the 
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strong energy condition. This means the Casimir force is attractive. But this is not always so in 
the cosmological context. Zel'dovich and Starobinskii [S] have indeed verified long ago that the 
Casimir energy of a scalar field could drive inflation in a flat universe with toroidal topology. 

The purpose of this letter is to exhibit a family of quantum scalar fields which originate a 
repulsive Casimir force in a closed universe, since they violate the strong energy condition. Inter- 
estingly, this family includes the simplest case one could consider: a minimally coupled massless 
scalar field with no potential. Our computation will be performed in the Einstein Static Universe 
(ESU), which can be faced as an approximation to a dynamical FRW model in sufficiently small 
time intervals, and avoids having to deal with the complexities of quantum field theory in a time 
dependent spacetime, like particle creation.^ 

2 Quantum scalar field with arbitrary coupling in the ESU 

We consider the ESU, which is a well known solution of the Einstein equations sourced by a perfect 
fluid with positive energy density (p > 0) and zero pressure (p = 0) together with a positive 
cosmological constant (A > 0): G^u + Ag^^i, = Tfj_^, with T^,y = pUf^u^. The metric is 

dslsu = -dt" + R'dns^ = -dt" + ^ ((4)2 + (4)2 + (4)2) . 

We have written the metric on the unit 3-sphere, dQgs, in terms of right forms on SU{2). In 
order to achieve a static solution, the cosmological constant and the energy density are related 
by A = p/2 = 1/B?. Another viewpoint is that the ESU is supported by a perfect fluid with 
p = -p/3 = -l/R^. 

It is well known that this universe is unstable against small radial perturbations as was first 
argued by de Sitter. The reason is that the energy density of the perfect fluid increases/decreases 
with decreasing/increasing radius, whereas the one of the cosmological constant is kept constant. 
Since the former gives an attractive contribution and the latter a repulsive one, any displacement 
from the original equilibrium position will grow, rendering such original position as unstable. But 
even without such classical perturbations this universe is unstable due to quantum mechanics. 
These are the instabilities we will focus on, in the spirit discussed at the end of the introduction. 

Let us consider a free (i.e with no potential) scalar field <I>, with mass p and with a coupling 
(not necessarily conformal) to the Ricci scalar of the background TZ = governed by 

(□ -^n)^ = fi^^ . 

Conformal coupling is obtained in four spacetime dimensions by taking the coefficient = 1/6 
(and the theory is then conformal if = 0), whereas minimal coupling corresponds to ^ = 0. The 
compactness of the spatial sections of the background guarantees a discrete mode spectrum which 
can be easily obtained using elementary group theory. We take the D'Alembertian in the form 

^Recent work on the Casimir effect on the ESU can be found in ^ and references therein. 
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where are the right vector fields dual to (t)j and is one of the two Casimirs in 50(4). Notice 
that the eigenfunctions of the Klein-Gordon operator (□ — ^TZ — fi^) may be taken in the form 

= e-^'^*p™^''"« , (1) 

where the index n represents all quantum numbers j, mi , m/j and '"^^ represents a Wigner 
D-function 0. Such function may be thought of as a spherical harmonic on the 3-sphere or as 
a matrix element of the rotation operator {j,mL\R{a, (3,^)\j,mii) , where \j,m) is the basis of 
a representation of SU{2) and (a, /3, 7) are Euler angles. It follows straightforwardly that the 
dispersion relation becomes 

with the degeneracy of each frequency being dj = {2j + 1)^, in agreement with the spectrum found 
in in]. Note that there are no unstable modes for ^ G W^, which includes minimal and conformal 
coupling. This is the range of couplings we will analyse in the following. 

Canonical quantisation of the scalar field can be performed unambiguously. One finds the mode 
expansion 



n ) 



with V = 27r2i?3 

being the volume of the constant t hypersurfaces and with the operators CLn,Q,n 
obeying the usual commutation relation [o„,a|j,] = 6nn'- 

The classical energy momentum tensor of the scalar field is more conveniently written in the 
natural tetrad basis = {dt, Ra}^/2, Ra'j^/2, Raj^/2}, 

2 

Tab = - ^k,$k^$ - ^5afe^' + CiGab " V^kfe + Qab^)^^ , (3) 

where we have denoted k^ = {d/dt, k^}. The conformal case = 1/6, /x = 0) has zero trace; quan- 
tum mechanically, however, the renormalised energy momentum tensor for a conformally coupled, 
free massless scalar field generically develops a trace anomaly, which can be written solely in terms 
of geometric quantities of the background jHUH]- In four spacetime dimensions such anomaly takes 
the form 

{T\)ren = J^^^ | C^afecd^'^'^' " J (iiabcd^?"'^' " ^RabR"" + 47^2)| + ^V'TZ . (4) 

The coefficient x is renormalisation scheme dependent. This trace is zero for the ESU: the Ricci 
scalar is constant, the Weyl tensor Cabcd vanishes since the geometry is conformally flat and the 
second Euler density also vanishes, since the Euler characteristic of any odd dimensional sphere is 
zero. 



3 Renormalisation 

Denoting the non-vanishing components of (T"^) as {T^ q) = —p and {T^ i) = {T"^ 2) — (^'^3) = 
we find the unrenormalised quantities 

1 1 +°° 2 2 

n=l n=l 
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with 



+ ^2i?2 +6^-1 



R 



n gN . 



(6) 



In [7j the particular case with ^ = 1/6 was studied. These infinite quantities were renormahsed by 
introducing a damping factor e~^" in the sums, subtracting the flat space contribution, and taking 
the parameter /? to zero. The result is 




n H 

2 8n 



41n^ + 1 + 47 



hoo +00 



n=Ar fc=3 



(7) 



1 



Pren 



127r2i?4 



Af-1 

E 



n=l 



1 + 

2 



3 (/^^)' 



3(^fi)^ 
8n 



n=Arfc=3 



Cfc- 



2fc 



n 



2fc-3 



1 , , (/.i?) 



+ 



+ 



121n^ +7 + 127 



120 24 32 

where 7 is Euler's constant, N is an integer such that > fiR and the explicit formulae for the 
coefficients bk, Ck are 

_ (-l)fe-i(2A:-3)!! _ (-l)^(2fc - 1)!! 

These expressions are not very enlightening; therefore we plot the renormahsed energy density and 
pressure, for flxed radius, in terms of the mass in figure H The point we would like to emphasise is 
that they are always positive. Therefore we conclude that a conformally coupled scalar field with 
arbitrary mass does not violate the strong energy condition and hence produces solely an attractive 
effect in the universe. 

At this point let us make a comment which also works as a consistency check. The first law of 
thermodynamics yields p = —dE/dV, where the total energy E = pV = p2iT^R^. We can easily 
check that this holds both for the infinite unrenormalised quantities © and the finite renormahsed 
ones (0 

1 d{R^po) 1 d{R^Pren) 



Po 



3i?2 dR ' ^"'^ 3i?2 dR 
Thus, we could have only computed the energy density and have deduced the pressure via the first 
law of thermodynamics, which obviously only states the conservation of energy T^'^;^ = 0. This is 
exactly what shall be done in the cases to follow. 

As a check on the result let us consider another renormalisation method, namely zeta 
function. We can write the regularised expression as 



where pq has dimension of mass and 



n 



47r2i?3-. („2 + q2)s/2 

n=l ^ ' 



(8) 



p'^R^ + ee - 1 , 
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Figure 1: Renormalised total vacuum energy 27r^pj.en and pressure 2'K'^pren for i? = 1, as functions 
of G [0,2], for conformally coupled scalar. Both quantities vanish asymptotically with increasing 
mass. 

which, for a? > 0, has the appropriate form to be written in terms of Epstein-Hurwitz zeta functions 

p{s) = 4^2^^ iCEH{s/2 - - a\EH{s/2,a^)) . (9) 

In order to obtain the renormalised energy density we should now consider the limit of this expres- 
sion when s ^ —1. 

The Epstein-Hurwitz zeta function is only defined by the infinite sum -|- a^)~'^ for 

Re{s) > 1/2. Its analytic continuation to other values of s, which converges for all s with the 
exception of an infinite number of poles is given by 

QEH{s,a^) = + V r(.) ^ + T{s) ^/'i^.-i/2(2vrnV^) , 

(10) 

where K^, are modified Bessel functions, and is required to be positive. It is simple to see that 
the infinite series converges for all s, since the modified Bessel functions have asymptotic behaviour 

Ku{z) ~ Y^e""^ , |z| ^ oo . 

The poles of CEH{s,a^) arise at the poles of T{s — 1/2), that is at s = 1/2 — n, n ^ Nq. 

Since a? must be positive this technique does not apply in the case of a massless minimally 
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coupled scalar field. Let us focus on the case with > 0. Applying H10() to @ we find 



' 2V I 4 r(f) + 
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(11) 

The point s = — 1 is exactly at one of the poles of the analytic continuation of the Epstein-Hurwitz 
zeta function. That is, the right hand side of the last formula is an infinite quantity, and it might 
seem that the zeta function method fails. With the correct interpretation this is not so. Let us 
consider separately two cases: 

• Conformal coupling = 1/6) and arbitrary mass: This case is obtained by substituting = 
fi^R^ in Hll|l . The infinite contribution, i.e the term with r(— 2) in the s = — 1 limit, is exactly 
the R independent term. Thus, despite its (infinite) contribution to the Casimir energy it 
will not contribute to the Casimir force. This argument would be enough to neglect it, and 
to suspect this is the flat space contribution. To confirm this is indeed the case compute the 
flat space result by taking the infinite radius limit of (ISJ: 

f+oD s+1 r-rfiziA) 

hm Pis) = ^ / dkk\e + M^)-^/2 = ^^Li22(^)3-. , 



R^+oc'^' ' An^ Jo ' ' 47r2 4 L (f ) 

This is exactly the infinite term obtained after zeta function regularisation. 

Nonconformal coupling: In the generic case, the divergent term cannot be identified with the flat 
space contribution, since it does depend on R. However, it turns out that the correct result 
is still obtained by simply dropping out the divergent terms, as will be shown below. Doing 
so, the final answer for the renormalised vacuum energy density of a scalar field with generic 
coupling is 



4.^ E ( ^2 i^-2(2W^^i2^ + 6e-l) 

n=l 
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An expression equivalent to H12|) was first derived, for the particular case of ^ = 1/6, in ^21) using 
a different renormalisation technique. It was argued therein that such renormalised energy momen- 
tum tensor could lead to a self consistent ESU solution of the semi-classical Einstein equations. We 
shall come back to this point in section 4. 

The advertised check on our calculation of the Casimir energy for a massive, conformally cou- 
pled, scalar in the ESU can be performed graphically; plotting the result (|12j) for ^ = 1/6, we have 
checked that it perfectly coincides with the previous one ((7|), plotted in figure ^ 

The zeta function method confirmed the damping factor method used originally in [7]. But is 
unfortunately inapplicable to the most general situation when a? <Q which includes the interesting 
case of a massless minimally coupled scalar field. We will now study this case by a variation of the 



6 



damping factor method. The essential ingredient wih be the Abel-Plana formula which we use in 
the form [121 



m=b 



^ ' 2 Jo exp{27rt) - 1 



(13) 



We regularise the vacuum energy density in the following way 

m 



^ +00 
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where is a positive function of its argument that grows sufficiently fast with n so as to make the 
sum converge for any (3 > 0. A convenient choice obeys 



it can be written explicitly as 
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The integration constant has been chosen such that 17 ^ as the frequency vanishes, which means 
that the damping factor is not altering the long wavelength modes. This seems to be the most 
physical choice, since these long wavelength modes are the ones responsible for the Casimir energy. 
We divide the situation with a stable spectrum into two cases which we analyse separately: 

• ^^i?^+6^ — 1 > (includes conformal coupling with arbitrary mass). Applying the Abel-Plana 
formula ([T3|l with 6 = to the regularised expression (fTl|) we find 
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The first integral became exactly the contribution of flat space (which justifies our choice of O). 
Indeed, from ((TH) 



R—f+<x 



+ 00 



lira p{p) = / (ixxVa;2 + /i2e-^^(^) = ^ / dOe"^^ , 



47r2 



+ 00 



(15) 



where 



0(x) = lim il(x) . 

R^oo 
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Figure 2: Renormalised total vacuum energy iTr'^pren and pressure 27r'^pren for i? = 1, as a function 
of ^ G [0, 1.5], for six different couplings ^ G [0, 1/6]. As ^ increases the colour of the line in the 
corresponding graph becomes darker. ^ = corresponds to the most negative curve in both graphs. 



The second integral, where a = \J /i^i?^ + 6^ — 1, converges in the /3 — > limit. Thus 



-co 



it+ia ^ ' y —it+ia 



exp{2'Kt) — 1 



dt . 



The square root representation in the complex plane that has been used, has a branch cut for 
t G [—a, a]. Splitting the t integral in + J^'^, one can check that only the second one contributes. 
The final result is 

1 r~ 'Vt^-(f^^fl^ (le) 



27r2i?4 yy^2^2+65_i exp{2TTt) - 1 

This expression was first obtained for the special case of conformal coupling ^ = 1/6 in lOj (see 
also j2]). Specialising for such coupling, the plot of this function exactly coincides with the ones 
in figure again checking the result. It is quite striking that expressions apparently as distinct as 
(O , (|T^ and (|16j) are actually different representations of the same function. 

For a more general coupling with //^i?^ + 6^ — 1 > 0, we can also check that the result (|16|) 
coincides with the one computed with the zeta function method ()12|). where we dropped out the 
divergent term, even though it could not be clearly interpreted as the fiat space contribution. Such 
term should seems to renormalise the gravitational action in the way described in ^1]. Since we 
are mostly interested in the renormalised energy- momentum tensor, which will be determined by 
the finite part of the result, we will not dwell any longer on this point. 

• — 1 < fi'^R'^ — 1 < (includes minimal coupling with zero mass). Applying the Abel-Plana 



8 



coupling 




-o.ooa 



Figure 3: Coefficient a in p^en = ajl^ and p^en = a/3i?^ for a massless scalar as a function of the 
coupling ^. 

formula H13|) with 6 = 1 to the regularised expression (|14() we find 

^^' '~Jn(±)4^^ 87r2i24 ^ ^ Att^R^ Jo exp{27rt) - 1 ' 

where 

The fiat space contribution is still (|T5|l . when //^ > 0. For the tachyonic case the fiat space 
contribution is more subtle, so we will restrict our analysis to — 1 < fi^R^ + 6^ — 1 < and jj? > 0. 
Subtracting the fiat space quantity and removing the regulator /? we find the renormalised energy 
density 
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(17) 

Despite the apparent unfriendly expression, this is a real quantity which can be plotted without 
great difficulty. Using the last formula and H16|) we have plotted several cases, with different ^'s, 
in figure [21 The most noticeable feature is that both the renormalised energy density and pressure 
may become negative for a range of values. Clearly this leads to violations of the strong energy 
condition, as will be discussed in more detail in the next section. Let us close by remarking that 
for zero mass, the renormalised energy density and pressure can be written, for arbitrary coupling, 
in the form 

a a 

Pren = j Pren = ^^^4 ) 
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Figure 4: Left: 27r^(pren + Sp^en) for six different couplings ^ G [0,1/6]. The colour scheme is 
the same as in figure When this quantity is negative the strong energy condition is violated. 
Right: 47r^(/>^g„ — p^en) fo^' C = 0; when this quantity is negative the velocity of sound exceeds that 
of light in the "quantum fluid" which implies (but is not necessary for) a violation of the dominant 
energy condition. 



where 



327r^ 



(60^ 



3V6^ - (6C 
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+ 



47r2 



/o(l + it)-/o(l-zt) 
exp{2TTt) — 1 



dt . 

(19) 

The quantity a is plotted as a function of ^ in figure |31 In particular, the special value of the 
coupling where a changes sign corresponds to a theory where the renormalised energy momentum 
tensor vanishes in the ESU. Its numerical value is — 0.05391. In figure |1] some energy conditions 
are displayed. 

4 Discussion 



The main purpose of this paper was to point out that the Casimir effect in the ESU becomes 
repulsive for a family of scalar fields with various couplings to the Ricci scalar and masses. This 
can lead to an inflationary era in the early universe, which generically seems to be too short to 
solve the usual big bang model problems. More interestingly it can lead to a cosmological bounce. 
To understand how, let us briefly consider the backreaction problem by taking the semiclassical 
Einstein equations, 

Analysing the simple massless case H18() . for q > 0, we conclude that the quantum fluid can support 
a self consistent Einstein Static Universe with some radius, a fact first noticed in J^]- Considering 
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the Priedmann equation and Raychaudhuri equations (reinserting Newton's constant) 

■o , 8-kG n ■■ AirG , 
R^ + k = ^pR\ R = —{p + 'ip)R, 

and taking the quantum fluid and a positive cosmological constant 

a , a 

P = A+^, p = -A + 



a self consistent solution is obtained with 



fc = +l, A = — , i?=^^_. (20) 

Clearly, this solution suffers from fine tuning, and is a universe of the order of the Planck size. 
For the case with a < 0, the fluid can, with the help of dust and a positive cosmological constant, 
produce an inflationary era followed by a decelerating phase and the present accelerating era. 
Indeed, taking 



one finds the generic behaviour displayed in figure El The universe has a minimal radius R = Rc 
where it bounces from a collapsing to an expanding epoch, thus avoiding the Big Bang singularity. 
For R > Rc there is a short accelerating phase, which is followed by a matter era and then 
cosmological constant domination. 

What is most interesting on having this Casimir energy induced bounce is that it does not 
rely on any particular form of a potential, and the scalar field needs not having classical effects 
whatsoever. It is its shear existence that originates the effect. Of course this is by no means a 
complete cosmological model of our universe. One obvious problem is how to include the radiation 
era in the picture, which has the same l/i?^ dependence for the energy density as the Casimir 
energy momentum tensor. Indeed, nucleosynthesis data constraints the influence that the quantum 
fluid could have over that period and therefore, since they vary equally with R, over any period 
where they both exist. A way around this problem could be to assume that radiation is only created 
after the quantum fluid domination epoch. Such mechanism could be similar to the usual reheating 
at the end of inflation. 

In any case, this simple model illustrates the point that the Casimir energy could both source 
an early inflationary epoch and avoid the Big Bang singularity in a closed universe. Note that 
despite having zero mass /i, the non-minimal coupling ^ works, in the Einstein static universe as an 
effective mass of order 1/ R^ . Thus, this massless case can have its classical dynamics frozen during 
inflation due to such effective mass. . 

It would certainly be interesting to further generalise this analysis to the case of a non-static 
FRW model where the dynamical Casimir effect takes place. 
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